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' Abstract 

CN 

, We study unitary random matrix ensembles in the critical regime where a new cut arises 

away from the original spectrum. We perform a double scaling limit where the size of the 
matrices tends to infinity, but in such a way that only a bounded number of eigenvalues is 
expected in the newborn cut. It turns out that limits of the eigenvalue correlation kernel are 
given by Hermite kernels corresponding to a finite size Gaussian Unitary Ensemble (GUE). 
When modifying the double scaling limit slightly, we observe a remarkable transition each 
time the new cut picks up an additional eigenvalue, leading to a limiting kernel interpolating 
between GUE-kernels for matrices of size k and size k + 1. We prove our results using the 
Riemann-Hilbert approach. 
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1 Introduction 

We consider unitary invariant random matrix ensembles on the set of Hermitian rax ra matrices, 
with a probability density of the form 

G\ . 

O. Z- 1 exp(-raTr^(M))dM, (1.1 

(N 



o 



where Z n is a normalization constant and dM is the usual flat Lebesgue measure on the Her- 
mitian matrices. We assume the confining potential V to be real analytic on 1R with enough 
growth at infinity, 

V(x) 

■ +oo as x —* ±oo. (1-2) 



log(x 2 + 1) 



Eigenvalues of a random matrix in the ensemble (jl.ip follow a determinantal point process 
generated by the following correlation kernel [22J, 



n-l 

K n (x,y) =e-5^)e-f^)^^ ) (x)^ ) (y), (1.3) 
given in terms of the orthonormal polynomials 

(n) / \ (n) k , M . n 

p K k J (x) = k\'x H , K y k ' > 0, 

with respect to the weight e~ nV on the real line. Using the Christoffel-Darboux formula, fjl .3|) 
can also be written in the following form, 

K n (x,y) = e-f^>e-f^> ^ (1 . 4 ) 



x - y 
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If we let the size n of the matrices grow to infinity, the limiting mean eigenvalue distribution 
of the ensemble exists and depends on V . In general it can be characterized as the equilibrium 
measure py (see e.g. [8j) minimizing the logarithmic energy in external field V, 

Iv{p) = jf log ^-^dp(x)dp{y) + J V(y)dp(y), (1.5) 

among all probability measures p on R. This minimization property is equivalent to the following 
Euler-Lagrange variational conditions [23]: there exists a constant £ G R such that 

2 J log \x — y\dpv(y) — V(x) = £, for x G supp py , (1-6) 

2 y log |jg — y\dpy(y) — V(x) < £, for x G R \ supp py. (1.7) 

It is known [10] that, for real analytic V, py has a density (py which can be written in the form 
<py{x) = Qy( x )i (1-8) 

where qy denotes the negative part of a real analytic function qy = qy — qy- Moreover qy(±x) 
is positive for large real x, from which it readily follows that supppy is a finite union of bounded 
intervals. The endpoints of the support are the zeros of qy with odd multiplicity. Generically 
the following conditions hold [19J: 

(a) the variational inequality (|1.7j) holds strictly for \ supply, 

(b) (py is positive in the interior of its support, or equivalently, qy has no zeros in the interior 
of supply, 

(c) (py behaves like a square root near the endpoints of its support, or equivalently, the zeros 
of qy with odd multiplicity are simple. 

In the critical cases where the above generic conditions do not hold, singular points occur. 
According to [TT], singular points are classified as follows. 

(i) Type I singular points or singular exterior points: these are isolated points outside supp py 
where equality in (jl.7p holds. Here qy vanishes at an order 4m — 2 for m > 1. 

(ii) Type II singular points or singular interior points: these are points in the interior of 
supp py where the density (py vanishes. Necessarily qy has a zero of multiplicity 4m at 
such a point. 

(iii) Type III singular points or singular edge points: these are endpoints of supp py where the 
density (py vanishes faster than a square root. Here the only possibilities are that qy has 
a zero of multiplicity 4m + 1 for m > 1 . 

When varying the potential V, the critical ensembles are the ones where a change in the 
number of intervals of supp py may possibly occur. Type I singular points correspond with the 
birth of a new cut away from the spectrum. Type II singular points indicate the closing of a 
gap in between two intervals of the support. Near type III singular points, a new interval can 
arise at the edge of the spectrum, or in other words, a gap can close simultaneously with one of 
the cuts. 
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Remark 1.1 It is important to note that not all multiple zeros of qy are singular points. For 
example, if supp/jy consists of one interval [a, b], it is known [8] that 

log\x-y\dp v (y)-V(x)-£ = -2 [ q 1 J\y)dy for x > b. (1.9) 



1/2 

Consequently, in order to have a type I singular (exterior) point x* > b, qy should necessarily 
change sign at some point in between b and x*. At this intermediate point, qy has a multiple 
zero although in general it is not a singular point. 

The local behavior of the eigenvalues of large random matrices near some reference point 
x* depends on the 'nature' of x* . Here 'nature' refers to the behavior of the limiting mean 
eigenvalue density near x* . The two regular kinds of points that occur, are points in the bulk of 
the spectrum (where ipy is positive) and points at the edge of the spectrum where cpy vanishes 
like a square root. The critical ensembles give lead, as described above, to three additional types 
of points that correspond to different local behavior of the eigenvalues. 

Local scaling limits of the eigenvalue correlation kernel (jl.3p turn out to be universal, which 
means that they depend on the nature of the reference point x* , but not on the confining 
potential V nor on the position of x* . In the bulk of the spectrum this leads one to the sine 
kernel [H EJ EU H21 123] , 

,. 1 v i * . u * , v \ smvr(u-tO 

J™. — 7Z^Z K n{x + — / ,v »S + — / ,vj - 



n— »oo 7Tipy(x*)n TTfy {x* )n ' 7T(py(x*)n 7r(u — v) 

Near a regular edge point, the limiting kernel is given in terms of Airy functions [9], 

1 Tr , * u * v , Ai (u)Ai'(v) - Aifv)Ai'(u) 
hm TnK n {x* + ^,x* + — ) = — ^ ^ ^ K -L, 

for some constant cy. 

As already mentioned, singular points indicate a transition where the number of intervals 
in the spectrum can change. These transitions can be observed when including a parameter in 
the confining potential, V = Vt- If the singular point corresponds to the value of t = t c , double 
scaling limits of the kernel, where we let n —* oo and at the same time we let t — > t c at an 
appropriate rate, lead to families of limiting kernels, depending on a parameter. Near singular 
interior points these double scaling limits are given by a kernel related to the Hastings-McLeod 
solution of the Painleve II equation, see [21 El [25] . In the singular edge case, the limiting kernels 
are related to a special solution of an equation in the Painleve I hierarchy [31 [7]. Singular exterior 
points where a new cut is born, have been studied in [16], but rigorous results about the limiting 
eigenvalue correlation kernel are not available in the literature yet. It is the aim of this paper to 
obtain rigorous asymptotics for the correlation kernel near a singular exterior point in the birth 
of a new cut. Unlike in the two other critical cases, there are no Painleve equations involved 
describing the local behavior of eigenvalues. In this case, the limit of the eigenvalue correlation 
kernel will be given by a kernel corresponding to a finite size Gaussian Unitary Ensemble (GUE). 
The size of the relevant GUE will depend on the precise choice of double scaling limit we take, 
or on the number of eigenvalues that are expected in the new cut. 

Similar transitions as the ones described above also occur in the study of the small disper- 
sion limit of the Korteweg-De Vries equation [18J, and in a more general context when studying 
Hamiltonian perturbations of hyperbolic systems [15[ [14] . In transitional regimes where alge- 
braic asymptotics for a solution of the KdV equation turn into elliptic asymptotics, phenomena 
are observed which are expected to correspond to the transitions for unitary random matrix 
ensembles, corresponding to the three types of singular points. Painleve asymptotics similar to 
those near singular interior points and singular edge points in random matrix ensembles have 
been verified numerically |T8], but not much is known about singular exterior points. This is 
yet another motivation to study the birth of a cut. 



3 



1.1 Statement of results 



The aim of this paper is to obtain a double scaling limit of the eigenvalue correlation kernel near 
singular exterior points. We deal with the case where qy has a double zero at x* , which is the 
lowest possible order of vanishing for a type I singular point. 

We consider a one-parameter family of potentials Vt = V/t, where V = V% is such that a 
singular exterior point x* is present, and such that supply = [a, b]. Furthermore we assume 
that there are, besides x* , no other singular points. For such a potential V, it was shown in 
pjj] that pt := pv t is supported on one interval [at,£>t] for t slightly less than 1, without singular 
points. It was also shown in this paper that for t slightly bigger than 1, a new cut is born 
near x*, so that supppt = [at,b t ] U [at,(3t\- As t \ 1, the cut disappears, so that at, fit — ► x*. 
The restrictions that supp py is supported on one interval, and that there are no other singular 
points besides x* , are technical rather than crucial. We expect that universality remains valid 
for potentials V that do not satisfy those assumptions. 

We work in a double scaling regime where we let the size n of the matrices tend to infinity, 
and at the same time we let t — ► 1 in such a way that 

„ , log n 
t - 1 < M— 2— 

n 

for some M > which can be arbitrary large. This double scaling limit implies that a bounded 
number of eigenvalues is expected in the vicinity of x* . 

Let us now formulate the main result of the present work. 

Theorem 1.2 Let V be real analytic satisfying condition M.ty) . and assume that x* is a type I 
singular point where qy has a double zero. Assume also that supply = [a, b] with b < x* , and 
that there are no other singular points besides x* . Let K n ^ be the eigenvalue correlation kernel 
tl.3\) for the potential Vt = V/t. We take a double scaling limit where n — > oo and t — ► 1 in such 
a way that \t — 1| < M-^|p. We define 

n f x * 1 

a:=2(t-l)- / - == ds, (1.10) 

; logn J b ^( s - a )(s-b) 

so that s remains bounded in the double scaling limit. Then, depending on the value of s, we 
have the following limits for the eigenvalue correlation kernel, 

V 1 K ( * U * V \ 

K GlJE (u,v;k) fork-\< s <k + \, k>l, 



for s < \ 



(1.11) 



2 ■ 



with c given by 



c=-Lq> v (x*) 1 /\ (1.12) 
and K GUE is given by 

K GU V v; k) = Jl e~^ as k > (L13) 

V 2 u — v 

where we write Hf. for the k-th degree normalized Hermite polynomial, with leading coefficient 

nk /2 2 

^ i/4 ^ , orthonormal with respect to the weight e~ x . 
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Remark 1.3 Obviously, the limits in ([l.lip do not hold not uniformly for s near a half positive 
integer. They do hold uniformly for s bounded and away from arbitrary small fixed neighbor- 
hoods of the half positive integers. 

Remark 1.4 K GUE (., .; k) is the eigenvalue correlation kernel for the k x k GUE, which is the 
random matrix ensemble for the potential V(x) = x 2 . After re-scaling, the eigenvalues in 
the newborn cut seem to behave asymptotically in the same way as the eigenvalues in a finite 
GUE. 

Remark 1.5 For s < 1/2, the limiting eigenvalue correlation kernel is trivial. This is not 
surprising since no eigenvalues are expected in the vicinity of x* for t < 1. A first eigenvalue 
near x* is only expected when s approaches 1/2. Each time we shift s with 1, an additional 
eigenvalue is expected in the new cut. 



Theorem 11.21 gives us the limiting eigenvalue correlation kernel in the case where s is not a 
half positive integer. It is natural to ask what happens when s is close to a half positive integer. 
It seems that, if s increases and passes a half positive integer, an additional eigenvalue is picked 
up by the new cut with high probability, which leads to a kernel corresponding to a GUE of 
larger size. Near the half integers, a remarkable transition takes place, involving a limiting kernel 
interpolating between a GUE kernel for matrices of size k and k + 1. 



Theorem 1.6 Under the same conditions as in Theorem there exist sequences X^ t such 
that the following asymptotics hold in the double scaling limit, 

K n i (x* + - " , x* + 



(cn)V2 n ' 1 \ {en) 1 / 2 ' {cn) l l 2 , 
' \-jK GVE (u,v;k) 

+A+,K GUE (u, V ;A: + l) + 0(^|), for k < s < k + 1, k > 0, 

Oin- 1 / 2 ), for s < 0. 



;i.i4) 



Furthermore the sequences X^ t are such that 



K,t + K,t — i> 

and 



K,t = 1 - K,t = 0(n^ 2+s - k ), as k < s < k + 1/2, (1.15) 

\- t = 1 - \+ t = 0(n 1/2+k ~ s ), as k + 1/2 < s < k + 1. (1.16) 

For simplicity in notation, we have written iT GUE (n, v;0) = in \1. 11$ . The expansion \1. 11$ 
holds uniformly for s bounded. 

Remark 1.7 Note that <|1.14[) is compatible with (II. lip . Indeed one observes that, for s away 
from a half positive integer, either of X~ t tends to 0, so that, in the limit, we are only left 
with one of the two kernels K GUE (n, v ; k) and K GUE (u, v; k + 1). This is exactly what is stated 
in Theorem II .21 Unfortunately, we are not able to give simple formulas for the sequences \ nt - 
We can give formulas for them, as we will do in Section [5j but we have been unable to reduce 
those formulas to simple expressions. 
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Remark 1.8 For s near a half positive integer, the limiting kernel changes abruptly. If we 
would put 



2 log n 

we expect that this would lead us to a limiting kernel of the form 

(1 - A)K GUE (u, v; k) + XK GVE (u, v;k + 1), 

with < A < 1. This is a kernel which corresponds to a determinantal point process as well, just 
like the GUE-kernels. The relevant point process is the one where we have k GUE-eigenvalues 
with probability 1 — A , and k + 1 GUE-eigenvalues with probability A. 

Example 1.9 A concrete example of a random matrix ensemble where a singular edge point is 
present, was given in [16J. For the potential 



1 ( T 4 e + e 3 ee — 2 2 



V ^ = rr^e U ~ + ~^r x + 2(e + ~ e)x ' ' (L17) 



with e > 2 and e such that fn( x — e)(x — e)y 'x 2 — 4dx = 0, there is a singular exterior point at 
x* = e, where qy has a double zero. For a potential of degree less then 4, singular points cannot 
occur. In order to construct an example for which there is a singular exterior point where gy 
has a zero of order 4m — 2, a potential V of degree at least 2m + 2 is needed. In those higher 
order cases, we expect that double scaling limits can be tuned in such a way that the limiting 
eigenvalue correlation kernel is no longer related to Hermite polynomials, but to polynomials 
orthogonal with respect to a weight of the form e~ p<yX \ where P can be any polynomial of degree 
2m. 



1.2 Outline for the rest of the paper 

In Section [21 we will construct equilibrium measures which correspond to the potential Vt- For 
technical reasons, we need modified measures compared to the usual equilibrium measures used 
in e.g. [8j HH [12], and also different from the modified measures used in |7]- In Section [3l 
we recall the Riemann-Hilbert (RH) problem for orthogonal polynomials introduced by Fokas, 
Its, and Kitaev [17]. We follow the ideas of the Deift/Zhou steepest descent method |13j in 
order to find asymptotics for the orthogonal polynomials. Here we follow similar lines as in 
[8j 111} [T2], with however two major differences. The first one is, as already mentioned, the use 
of modified equilibrium measures, and the second one is the construction of a local parametrix 
near the singular point x* . For this construction, we will need, in Section 01 a model RH problem 
built out of Hermite polynomials. The local parametrix will enable us to find asymptotics for 
the orthogonal polynomials near x* . In Section \5\ we will use those asymptotics to obtain 
asymptotics for the eigenvalue correlation kernel and to prove Theorem 11.21 and Theorem 11.61 



2 Equilibrium measures 

It is a well-known fact that a ^-function related to an equilibrium measures plays a crucial role in 
the Deift/Zhou steepest descent analysis. As already mentioned in the introduction, the limiting 
mean eigenvalue distribution pt is an equilibrium measure in external field Vt, and would hence 
be an obvious candidate to built out the (^-function. However this would not be a convenient 
choice because the endpoints at and /3t of the new cut vary with t and both tend to x* as t \ 1. 
This would create various technical difficulties for the construction of a local parametrix near 
the critical point x* , which will be the most crucial issue in our RH analysis. 
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To prevent the presence of endpoints near x* varying with t, we will construct for t > 1 a 
modified equilibrium measure pt which we force to have its support away from x*. The portion 
of mass of p t near x* can however not be ignored, and for this purpose we add to p t a point mass 
centered at some point close to x* . The rough idea is that we can 'approximate' the limiting 
mean eigenvalue density p t by a probability measure of the form p nit = p n j + m n ^8 x * , where 

• p nt t is the positive equilibrium measure in external field Vt with mass 1 — m n t, with m n t 
given by 



m n t = max 



2 t-l_ rx* — i d as t > 1, 

logn Jb ^/( s _ a )( s _b) ' ^ 1) 

as t < 1. 



b n,t 



is the Dirac distribution centered with mass 1 at x* t , 
is a point near x* which we will determine below. 



To be precise in defining p ni t, it is the equilibrium measure minimizing the logarithmic energy 
Iv t (p), defined by (jl.5j> . among all positive measures p satisfying the following two conditions: 

• supp p C R \ [x* — e, x* + e] for some sufficiently small fixed e > 0, 

= 1 — m n .t, with m n ^t defined by (|2.1j) . 



The equilibrium measure p n t does not depend on the choice of e, at least not if e is sufficiently 
small so that [x* — e, x* + e] does not intersect with [a, b], and if t is sufficiently close to 1. Note 
also that, for t < 1, our new equilibrium measure is exactly equal to pt- The measure p, n t is 
thus independent of n for t < 1. 

Since there are no singular points in external field V any longer after the exclusion of [x* — 
e,x* + e], it was shown in [19] that for t sufficiently close to 1, supp/i n; ( : — [fnii^t] consists 
of one single interval. It is a standard fact [23] that p n j satisfies the following variational 
conditions, with Vt = V/t, for t sufficiently close to 1, 

2 J log \x - y\dp ntt (y) - V t (x) = t n j, for x £ [a' n<t , tf n>t ], (2.2) 

2 / log \x - y\dp nt t(y) - V t (x) < l n , t , for x G M \ {[a' n ^ V t ] U [x* -e,x* + e}). (2.3) 



Those variational conditions will be crucial throughout the following sections. 

We can directly apply the results obtained in [10] to conclude that the density V'n.t of p Ut t 
can be written in the form 

^t(x) = -jQ- t (x), (2.4) 
where Q n j is a real analytic function with negative part Q^t- Also it follows from [10] that 

feW = fmv '/mi!i^ w (2 . 5) 



2t ) t _ 

Using this identity, weak* convergence of p, n ^ to p\ = p n ^\ is enough to conclude that Q n ,t(z) — ► 
gy(z) uniformly on compact sets in a neighborhood of the real line, as t — > 1 and n — > oo. This 
means that Q n j has simple zeros a' nt and tending to a and 6ast— ► 1, n — >oo. However 
we can do more. It follows from a result by Buyarov and Rahmanov [4] that 

Mn,t-^n,l = 0(t-l), (2.6) 
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Using (|2.5|) we now have that Q n ,t( x ) ~ Qv( x ) = C(t — 1) uniformly on compact sets, which 
implies that a! n t = a + 0(t - 1) and b' nt = b + 0(t -1). 

It is convenient to rewrite (|2.4[) in the following way, 



where h n t is real analytic on R and h n ^ = h n \ + 0(t — 1) uniformly on compact sets as t — > 1, 
n — > oo. Now it follows from formula (jl.8p and the fact that qy has a double zero at x*, that 
h n \ has a simple zero at x*. It then follows that h n ^ must have a zero x* 4 = x* + 0(i — 1) as 
t — > 1, n — > oo. This zero x* 4 is the point where we center the Dirac measure approximating 
the portion of the limiting mean eigenvalue distribution in the new cut. Note that the point x* t 
is a double zero of Q n ,t, but it is not a singular point, cf. Remark II .li 

So far, we did not give any arguments why we need to choose m n ^ and x* t in the way we 
did. Only at the very end of the RH analysis, in the construction of a local parametrix near x* 
in Section [U it will become clear that only these choices for m n ^ and x* t do the job. Except 
for the construction of the local parametrix, the RH analysis would also work for other values 
for m n ^ and x* jt . 

3 Riemann-Hilbert analysis 

The starting point of our analysis is the RH problem for orthogonal polynomials introduced by 
Fokas, Its, and Kitaev [17] . We follow the approach of [8j [Tl], [12] , where the Deift/Zhou steepest 
descent method |13| was used in order to find large n asymptotics for the solution of this RH 
problem. The main idea is to apply a series of transformations to the RH problem in order 
to find, at the end, a RH problem which can be solved approximately for large re. The crucial 
new feature in our analysis is, besides the modification of the equilibrium measure described 
in the previous section, the construction of a local parametrix near the singular exterior point 
x* in a double scaling limit. This will be done in Section [J] using the RH problem for Hermite 
polynomials. The remaining part of the Deift/Zhou steepest descent analysis follows similar lines 
as the ones developed in [51 [Til 12] and later also applied to double scaling limits in [51 El Ej • 

For technical reasons, we assume that the confining potential V is such that the support of 
the limiting mean eigenvalue density consists of one single interval, supppy = [a, b}. We also 
assume that the singular point lies at the right side of the spectrum, x* > b, which we can do 
without loss of generality because of the possibility to consider the potential V(—x) instead of 
V(x). Another technical restriction is that we assume the absence of any other singular point 
besides x* . 

3.1 RH problem for orthogonal polynomials 

For each n and t, we consider the following RH problem. We seek for a 2 x 2 matrix- valued 
function Y(z) = Y(z;n,t) satisfying the following conditions. 

RH problem for Y: 

(a) Y : C \ R -» C 2x2 is analytic, 

(b) Y has continuous boundary values Y±(x) for x G R, where Y + (x) and Y-(x) denote the 
limiting values when approaching x from above and below, and 



lf>n,t( x ) = ~\[Wn,t - x )( x ~ a 'n,t) h nA x )X[a' n t ,b' n j(x) 



(2.7) 




e -nVt(x) 
1 




(3.1) 
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(c) Y has the following asymptotic behavior at infinity, 

1 



Y(z) = (l + 0(z^)) 



Z~ 



as z — ► oo. 



(3.2) 



This RH problem has a unique solution which is given in terms of the orthonormal polyno- 
mials pk = pi with respect to the weight e~ nVt on M, see |17j . 



Y(z) 



( 



X 



'"n 

2tH 



K n Pn( z ) 

>,*) 



Pn{u)e 



-nVt(u) 



du 



\ 



u — z 

p n -l(u)e~ nvt ^ 

u — z 



du 



for z G C \ M. (3.3) 



Here we have written = K k > for the leading coefficient of p^ . 

It is straightforward to check that the eigenvalue correlation kernel K n j, given by (|1.4j) . can 
be expressed in terms of Y. Using the fact that dety = 1 (this follows from the RH conditions 
for Y using a standard complex analysis argument), one obtains 



K n , t (x,y) = e -f^) e -fv t (,)__I__ (0 !) Y ^(y)Y ± ( X ) Q . 



(3.4) 



From (|3.4p . it is clear that all the information needed to prove Theorem 11.21 is contained in 
the RH solution Y. We now need to find sufficiently accurate asymptotics for Y in the double 
scaling limit. 



3.2 First transformation Y i— ► T 

In this section we perform a transformation Y t— ► T of the RH problem which normalizes the 
behavior of the RH solution at infinity in such a way that T(z) — ► I as z — > cxd. Besides that, the 
map will modify the jumps in a convenient way. A crucial role in the transformation is played 
by the so-called ^-function, which is related to the equilibrium measures constructed in Section 

El 

The standard way to define the ^-function would be to put it equal to J log(z — s)dpt(x), 
where pt is the limiting eigenvalue distribution of the random matrix ensemble. As already 
noted in the previous section, this would involve several technical difficulties. We have already 
anticipated to these difficulties by replacing the limiting mean eigenvalue distribution pt with a 
modified measure p Ut t = p n ,t + fnn,t$x* n t , where we approximated the portion of mass in the new 
cut by a Dirac measure centered at x* t . The new equilibrium measure leads us to the following 
^-function, 

9(z) = 9n,t( Z ) = J lo g( Z ~ y) d Pn,t(y) = J 1 °&( Z ~ V) d ^n,t{y) + m„,t log(z - Z* )t ), (3.5) 

where we take log z analytic in C \ (— oo, 0] with — ir < Im log z < ir. The variational conditions 
(12. 2|) and (|2.3I) can now be translated to properties for the g-functions. For t sufficiently close 
to 1, we have that, with £ = £ n j, 

9+{x) +9-{x) - V t {x) - £ = 2m njt log|x* jt - x\, for x e [a' n)t , b' n>t ], (3.6) 

g+(x) +g-{x) - V t (x) -£ - 2m„ )t log|a;* jt - x\ < 0, 

for x G R \ (K )t) b' n j] U [x* - e, x* + e]). (3.7) 
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Note also that 

r+oo 

9+{ x ) ~ 9-(x) = 2m / dfi n ,t(y), 

J x 

which means in particular that 



for x G 



(3. 



g+(x) - g-{x) = < 



2m\ for x < a' n t , 

27rim n j, for b' n t < x < x* t , 
0, for x > x* nt . 



(3.9) 



The above properties of the ^-function are, together with the fact that e n9<yZ ^ = z n (l+0(l/z)) 
as z — > oo, crucial to transform the RH problem. We define T as follows, 



T(z) = e-¥^Y{z)e- n9{z)a *e^ a \ for z G C \ M, 



(3.10) 



where 03 is the third Pauli matrix, o 3 = {\ \ ' With v = r 11 ^! 5 ) 0} = nm n,t, it is straight- 
forward to check using (|3.6p . (|3.7|) . and (|3.9p . that T satisfies the following RH problem. 



RH problem for T: 

(a) T : C \ R C 2x2 is analytic, 

(b) 7+(x) = T_(x)vt(x) for x£l, with 

f/l e »C0+(aO+0-(*)-Vt(*)-*)' 

,0 1 



i>t(x) 



-n(g+(x)-g- (x)) 




x — X 



* \2u 
n,t\ 



-2-kiv p n(g+(x)+g-(x)—Vt(x)- 







for x e (-00, a^ >t ) U (x* >t , +00), 
forxG K jt ,6' nit ], (3.11) 
for x E (6' t ,x* t ), 



(c) T(z) = I + C(l/z), as z 00, 

(d) T{z){z — x* J 1 ^ 3 is bounded near x*. 

It is necessary to add condition (d), controlling the behavior of T near x* , in order to have 
unique RH solution. A similar condition was not stated in the RH problem for Y, where we 
assumed continuous boundary values and thus a bounded RH solution near x* . In exception of 
x* t , T has continuous boundary values on R as well. 

Using (|3.4p and (|3.10p . we find the following identity for the eigenvalue correlation kernel in 
terms of the new RH solution T, 



K ntt (x,y)=e-n Vt We 



2^7)(° l >WV<*)(J).w.m 
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3.3 Second transformation T i— »■ S 

Before going on with the next transformation of the RH problem, we first write the jump matrices 
vt in a slightly different and more convenient way. Let us define a function <j> = <j> n t, analytic 
in C\ [-00, V t ], by 



<j>{z) = f ™ &l 2 t (s)ds = f n '\s - a' n ^l\s - b' n ^ 2 h n , t (s)ds, (3.13) 

J z J z 

with Q n j as in (|2.4p and h n> t as in (|2.7p , and with a suitable branch of the square root such that 

Ke(j)±(x) < 0, as x < a' nt , 

±Im ± (x) > 0, as x G [<4 )t , , 

<^(x) < 0, as x G (^j, +00) \ [x* — 5, x* + 8]. 

Using the definition of (f>, (|2.7p . and ()3.9j) . we observe that 

9+{x) - g-{x) = ±24>±(x) + 2mm n!t , for x G [a' n<t , b' n>t }. (3.14) 
On the other hand, by (|3.6p . we have that 

5±0*0 - 5t( x ) = 2 9±{x) - V t (x)-£ - 2m njf log|x* jt - x\, for x G K t ,6^J. 

Combining the two above equations and using the identity theorem, we obtain the following 
identity, 



2<t>{z) + 2m nii log(2 - < it ) = 2 5 (z) - F t (z) - £, for z G C \ (-00, a?* t )- 



(3.15) 



We can now rewrite the jump matrix vt given by (|3.1ip . using f)3. 14|) and (|3.15p . For t > 1 
we have the following jump matrix, 



' / 1 e 2 ™|x — x* | 2l/ e 2n ^+( x ) N 



ut(x) = < 



and for t < 1, 



■yr(x) = < 



q e -2n</>_(x) e 27rw 
2i/ e 2n^)(a;)\ 



e 2mv \x - .r 



71. t 

2iviv 



1 |x-x; it | 2iy e 2n ^) N 



for x G (—00, a' n)t ), 
for x G K it ,?4 >t ], 
for x G (^ jt ,x* t ), 
for x G (x* t , +00), 



' e -2n</> + (x) J 

e -2n<p-{x) 

\ e 2n4>{x)\ 

, • 



for x G (a' t ,b' t ), 



for x G (-00, a' t ) U (6! ' t , +00) 



(3.16) 



(3.17) 



This way of writing down the jump matrix is convenient because it is now clear that we can 
factorize the jump matrix on [a' nt ,b' n t ]. Indeed, for x G [a' nt ,b' nt ], the jump matrix can be 
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written as follows, 



v T {z) 



x — X. 



* \-1v 
n,t\ 





g— 2n(fr—(x) g2iriv ^ 



X — X 



* 1 2u 



\X — X 



* \ — 2v 
n,t I 







x — x. 



* I — 2v p — 2n(f>+(x) „— 2mu i 

. (3.18) 



Because we can extend the first and the last factor analytically to the lower resp. the upper half 
plane, this factorization allows us to deform the jump contour in such a way that the different 
factors of the jump matrix lie on different curves. This is very convenient because it allows us 
to transform jumps that are oscillatory with n on the real line to jumps that are exponentially 
decaying on some contour in the complex plane. We deform the jump contour to a lens-shaped 
contour Eg as shown in Figure [TJ 



We define 



S{z) 



{ T(z) 
T(z) 

T(z) 



0' 

— 2-kiv 1 



outside the lens-shaped region, 
in the upper parts of the lens, 



{X* nt - z )-2v e -2n<t>(z) e 2iriv j 



in the lower parts of the lens. 



(3.19) 



Due to the opening of the lens, the analytic continuations of the three factors of the jump matrix 
vt now live on different contours. We obtain the following RH problem for S. 

RH problem for S: 

(a) S : C \ S 5 -»■ C 2x2 is analytic, 

(b) S+(x) = S-(x)vs(x) for with 



' / 1 e 27Til/ \x — x* t \ 2l/ e 2n ^ + ^ x ^ 



v s {x) = < 







1 



1 

{z-x* nt )- 2v e- 2n ^ 1 r 





\x — X 

-2 





\x — X 



* 1 2u 
n,t I 



n.t I 



„— 2-kiv |„ \2v „2ncf>(x) 

„2niu 



1 \x-x* nt \ 2u e 2n ^ 



for x G (-00,0^), 
for x £ SsnC*, 
for x £ [a' n)t ,b' n)t ], 
for x G {b' n ^x* n t ), 
for x G (x* n , t ,+oo), 



(3.20) 



(c) S"(z) = J + 0(1/^), as z -> 00, 

(d) S(z)(z — x* nt ) uaa is bounded near x* 
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Figure 1: The lens-shaped contour Eg 



Using the fact that S(z) = T(z) for z outside the lens-shaped region, and (|3.15p . we find 
from (|3.12p that, for x, y in some sufficiently small fixed neighborhood of x* , 



Kn,t^y) = - n 2 m{ n x _ y) (0 l)S^(y)S + ( X )^ Q j. (3.21) 

3.4 Construction of parametrices 

If we take z away from arbitrary small but fixed neighborhoods surrounding a, b, and x* , it 
can be checked as e.g. in [8j using the Cauchy-Riemann conditions that, for n — * oo, the jump 
matrix vs(z) converges exponentially fast to a jump matrix v^°°\z) given by 



v^(x) 



\ x x *n,t\ 2U 

-\x-x* n ^ J 



for x £ [a' n<t ,b' n j], 

for x G (b' n:t ,x* n t ), 
elsewhere. 



(3.22) 



In other words, the jump matrices on the lips on the lens and on R\ [a— e, x* +e] are exponentially 
close to the identity matrix. As t > 1, on (6 + e, x* — e), in the gap in between the two intervals 
of the spectrum, the off-diagonal entry is also exponentially small. 

If we ignore, for a moment, small neighborhoods U a , Ut>, and U x * surrounding a, b, and x*, 
the RH problem for S is reduced, up to exponentially small jumps, to a RH problem which 
we call the RH problem for the outside parametrix, referring to the region away from the 
local neighborhoods of the special points. Besides the outside parametrix, we will need local 
parametrices near a, b, and x* points in order to obtain uniform asymptotics for S. The outside 
parametrix will determine the asymptotics for S away from the special points, while the local 
parametrices will contribute to the local behavior of 5. Of particular importance for us is the 
construction of the local parametrix near x* in Section HJ which will in the end describe the 
local behavior of the eigenvalues in the new cut. 

3.4.1 Outside parametrix 

Ignoring the exponentially small jumps and small neighborhoods of the special points, our RH 
problem reduces to the following. 

RH problem for p(°°) 
(a) P(°°) : C \ [a' t ,x* t ] -► C 2x2 is analytic, 
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(b) p(°°) satisfies the following jump conditions, 

PtH*)=P^\x)(_ lx _ <tr2u \ x ~<f U y as *G (3-23) 
PtH*) = Pt\x)e- 2 ™°s, as x G (^ t ,< t ), (3.24) 

(c) P(°°)(ar) = 7 + 0(l/z), as z -> oo. 

It is clear from the jump conditions that p(°°) = P^^ will have some singular behavior near 
a 'n ti Kit) an d x *n,t- Since we did not specify the required behavior near those points, the solution 
of this RH problem is not unique. However there is only one solution which is compatible 
with the local parametrices that we will construct afterwards. We now construct this solution 
explicitly. 

lit < 1, we have that v = so that p(°°) is analytic in C \ [a' nt , b' n t ]. One checks directly as 
e.g. in [8] that 

P {oo) (z):=[ m+ % t) -i ) , *eC\[< t ,5U (3.25) 

V 2i 2 / 

with 

(z-V \ 1/4 

0(z) = J3 n>t (z) = [ ^ , zeC\[a'b']) (3-26) 



z — a. 



n,t 



is a solution of the RH problem for the outside parametrix if t < 1. 

For i > 1, the situation is slightly more complicated because of the additional jump on 
(b' n t , x* t ) . This jump can be created by introducing an auxiliary scalar function D, which is 
analytic in C \ [a' n t , x* n t ] and has the following jumps, 

D + (x)D_(x) = \x-x* nt \ 2u , for x G [a' nt ,b' n ^, 

D + (x)D_(x)- 1 = e 2 ™, for x G [6^,< >t ]. 

In addition the limit 

Doo = lim D(z) G R, (3.28) 

z— too 

should exist, so that we can define p(°°) as follows, 

p(°°)(z) = D^p(°°)(z)Z)(z)- ,T3 . (3.29) 
Using (|3.27p and the fact that 

pf ) (a .) = p0») (a .) ^ J) , for x G « t> 6^), 

one verifies that indeed p(°°) is a solution of the RH problem for the outside parametrix. 
We will now construct the function D. Let us first consider the function 



$(z) = z+y/{z-l){z + l) for z G C\ [-1,1], (3.30) 
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which is the conformal mapping from C \ [— 1, 1] to the exterior of the unit disk. This function 
has the convenient property that $ + (x)<I>_ (x) = 1 for x £ [—1,1]. If we let F map [a' n t ,b' nt ] to 

[-1,1], 

u n,t u n,t u n,t u n,t 

we have that $ + (F(x))$-(F(x)) = 1 on [a' n ^b' n ,t]- 
We also need the function 



^(z-<,t)(z-b' n ,t) fK,t \og{x* nt - X ) 



dx 



z — X 



G(z) = exp V / 

\ * J << yl(*-<,tW n ,t-x) 

which is analytic in C \ [a' n t , b' n t ], and using a residue argument we find that 

G + (x)G-(x) = x* n t - x, for x £ [a' n ^b' n>t \. 

Now let us write v = k + A, with k £ N U {0} and |A| < 1/2, so that |A| is the distance 
from v to its nearest nonnegative integer. Define D for z £ C \ [a' n t,b' nt ] in the following way, 

D(z) = (z- x* nit ) A <S>(F(z)y A G(z) k . (3.32) 

Since the branch cuts of the first two factors cancel out against each other on (— oo, a! nt ), D is 
analytic in C \ [a' n ti x nt\i an d one also checks that D has a limit as z — » oo. Using the jump 
properties of &(F), G, and (z — x* n t ) A on {a' nt ,b' nt ) and on (b' nt ,x* nt ), one verifies that (13.270 
is satisfied, so that p(°°) solves the RH problem for the outside parametrix. 

As mentioned before, has some singular behavior near a' nt , b' nt , and x* nt . We have 

that 



n,ti 



P^\z) = 0(\z-a' n ,\- 1/4 ), 
P^\z) = 0{\z-b' ntt \- 1/A )i as^l^ 
PN( Z ) = 0(\z-x* nt \-\ A \), asz^i 



n,t ' 



This is the 'least singular' behavior that p(°°) can possibly have if it satisfies the prescribed 
jump conditions. 

3.4.2 Local parametrices near a and b 

The construction of the local parametrices near a and b uses a model RH problem with a solution 
built out of the Airy function and its derivative. This construction can be done in exactly the 
same way as in [8] for v = and as in |21| for v > 0. The precise construction of the parametrix 
is not important for us. At this point, it is enough to know that local parametrices in sufficiently 
small but fixed neighborhoods U a and £/& of a and b (note that a' n t and b' n t are included in those 
neighborhoods for t close to 1 and n large) exist in such a way that 

(a) P : U a U U b \ S 5 -» C 2x2 is analytic, 

(b) for z e £ 5 n (J7 a U £/" 6 ), we have P + (z) = P-(z)vs(z), 

(c) for z £ dU a U dUb, if we let n — > 00 and t — > 1, we have 

P^P^)^)- 1 = / + 0(n- 1 ). (3.33) 

This matching of P with p(°°) can only be obtained because the behavior of p(°°) near the 
endpoints a' nt and is 'not too bad'. If we would have chosen a different outside parametrix 
with different behavior near the endpoints, this would not be the case. 
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Figure 2: The contour £# after the third and final transformation. 



3.4.3 Local parametrix near x* 

The crucial part of the RH analysis consists of constructing a local parametrix in a sufficiently 
small but fixed neighborhood U x * of the singular exterior point x* . Our aim is to find a function 
P satisfying the following conditions. 

RH problem for P 

(a) P : U x * \ £ 5 -> C 2x2 is analytic, 

(b) P+(z) = P-{z)v s {z) for z E U x * n £<?, 

(c) If we take the double scaling limit where we let n — > 00 and at the same time we let t — ► 1 
in such a way that |t — 1| < M^^, and if we then put 

n r x * 1 

a = 2(t-l)- / — == da, (3.34) 

we have that 

I + C(n-2+l A l) for z E d£4», if s > 0, 
P(s) = p(°°)(z) x { t (3.35) 

I + C(n~2) for z € SCTa-., if s < 0, 

where once again we have written v = k + A, with A; E N U {0} and |A| < 1/2. 

We postpone the construction of P to Section |H We assume for now the existence of the 
parametrix P satisfying the above conditions, and proceed with the remaining part of the RH 
analysis. It is important to note that the construction of the local parametrix near x* will only 
work for our particular choices of m n t and x^ t in the construction of the equilibrium measure. It 
is remarkable that this is the only part of the RH analysis that would fail for arbitrary bounded 
v or for an arbitrary sequence x^ t tending to x* . We note already that, in (|3.35p . the 'matching' 
of P with p(°°) breaks down if A = ±1/2, or when v is a half integer. 

3.5 Final transformation of the RH problem 

We now define the function R as follows, 



{ S(z)P^(z)-\ for z outside the disks. 

Here P is the parametrix satisfying the RH problem posed in Section 13.4.31 inside U x * , and 
P is equal to the Airy parametrices described in Section 13.4.21 inside U a and [/&. Inside the 
disks, the parametrices are constructed in such a way that they have exactly the same jumps 
as S has, and this implies that R is analytic inside U x *, U a , and £/&. R is also analytic on 
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\a! n i, b' n J \ (U x * U U a U Uf)) since S and P(°°) have the same jump here. We can conclude that P 
is analytic outside a contour as shown in Figure [21 Outside the disks, the jump matrices vr 
converge exponentially fast to the jump for p(°°). At the boundaries of the disks, the jump 
matrices converge as well because of the matching of the local parametrices with the outside 
parametrix. One verifies that, with the contour orientated as indicated in Figure [21 R satisfies 
the following RH problem. 



RH problem for R 

(a) R : C \ E R -> C 2x2 is analytic, 

(b) R+{z) = R-(z)vr(z) for z G Hr, with vr given by 

l+0(e~ cn ), forzGS R nS 5 , 
P(z)P^°°\z)- 1 , for z G dU x * U dU a U 9?7 6 . 



(3.37) 



(c) R(z) = 1 + 0{z~ l ) as z -> oo. 

In the double scaling limit where we let n — ► oo and at the same time we let t — * in such a 
way that \t — 1| < M^^ 1 , we recall from (13.330 and (I3.35P that 



P(z)P {oo) (z)- 1 



' / + C(n- 1 ), for z e dU a UdU b , 

I + C(n- 1 / 2 +l A l), for z G 9f7 x * , if s > 0, 
, / + C(n- 1 /2) ) for z G 9{7 X * , if s < 0. 



(3.38) 



Using (|3.37p . we see that the jump matrix vr has the following asymptotics in the double scaling 
limit. 



vr{z) 



I + 0(e~ Cn ) 



for z G D Us, 
for 2 G dU a U aC/ fc 



J + O^- 1 / 2 ^^), for z G SC^., if s > 0, 
/ + 0(n~ l l 2 ) for z G 917a,., if a < 0. 



(3.39) 



Here C is an unimportant positive constant. 

From the uniform convergence of the jump matrix to the identity matrix (except when 
A = ±1/2), it now follows as e.g. in 0Q2] that in the double scaling limit, also the RH solution 
R is uniformly close to the identity matrix, 



R(z) 



7 + 0(n- 1 /2+|A| )i if a > o, 
7 + C(n~ 1 /2) if s <0. 



(3.40) 



In fact, one can even obtain a full asymptotic expansion in negative powers of n, but for our 
purposes, (|3.40p is already sufficient. These asymptotics hold moreover uniformly for z G C\£r, 
and uniformly in s as well, although for s near a half positive integer, (|3.40p only says that R is 
uniformly bounded, without convergence to the identity matrix. 
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4 Construction of the local parametrix near x* 



In this section, we construct the parametrix P solving the RH problem posed in Section 13.4.31 
explicitly. The construction is different in the case where t < 1 compared to the case where 
t > 1. For t < 1, the support of \x n ^ consists of one interval, and the parametrix can immediately 
be written down in terms of elementary functions. For t > 1, we will use a model RH problem 
built out of Hermite polynomials. This is the crucial new feature in the RH analysis that makes 
it possible to find double scaling asymptotics near the singular exterior point. 



4.1 Construction of the parametrix if t < 1 

Recall that we are in the one-interval case with v = for t < 1, and that the equilibrium measure 
does depend on t but not on n in this case. For convenience, we drop the subscript n to write 
the functions, related to the equilibrium measure, for which it is clear that they do not depend 
on n if t < 1. 

We can construct the parametrix P directly as follows, 

c 2n<A t (a:*) 



( e *n<l> t (x t ) r. e ~" J \ 



(4.1) 



Here we should still define ft as a conformal mapping in U x * , positive for x > x* n t , in such a way 
that P has the appropriate jump condition. Using Cauchy's theorem and the fact that is 
analytic in U x * , one verifies that 

/1 2n4> t (x* t ) p -nft(x) 2 \ 

P+(x) = P-{x) ( J * j , for x E Ua* n M. 

As explained in Section 13.4.31 we want to construct P in such a way that 

/1 2n0 t (x)\ 

P+(a:) = P_(x)«5(x) = P_(a:)fJ x , for x G 17*. n R, 

which is now equivalent to the following condition on ft, 

ft(z) 2 = 2Mx*t) ~ 20t(«)- (4-2) 

Since Xt is defined as the simple zero of h t near x*, we have by (|3.13j) that (f>t(z) — 4>t{x* t ) has a 
double zero at x%, and consequently ft, defined by condition (|4.2p . is a conformal mapping near 
x\. Furthermore we have that 

ft(x* t ) = 0, fl(x* t ) = ±- A Q'l(xl)^ > 0. (4.3) 



It remains tho show that the matching condition (c) in the RH problem for P is valid as 
well. For this purpose, we note that 4>t( x t) < if t < 1. Indeed this is the case because we know 
by flSIgD that 

Mx*t) = 2gt(x* t )-Vt(x* t )-i t . 

The right hand side is negative for t < 1, since the opposite would violate the variational 
inequality (|1.7j) for the equilibrium measure fit = Pt- By (|4.1I) we now find easily that 

P(z) = P {oo \z)(I + 0{n- 1 ' 2 )), asn^oo, t < 1. (4.4) 

This shows that P satisfies the required RH conditions, stated in Section [3.4.31 
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4.2 Construction of the parametrix if t > 1 

Because we do not have a variational inequality for [x n t near x* any longer if t > 1, the con- 
struction of the parametrix, as done in the previous section, fails for t > 1. In this case we have 
that v = s > 0. 

4.2.1 Model RH problem for * 

In order to construct the parametrix P, we need the function = ^(C; k) defined for k £ NU{0} 
by 



*(C;ft) 



2 (fc-l)/2 



Tr^Vfc! /" H k (u)e- u2 



7T 1 /4 > /( fc _l)! 



2 • 2 fc / 2 7Ti 
2 (fe-l)/2 



(if/ 



\ 



7T 1 /4 N /(fc_ 1)! 



u-C 

g t -i(Me-"' 
u-C 



for c e c \ : 



(4.5) 



where denotes the degree k Hermite polynomial, orthonormal with respect to the weight 
e~ x on R. It is a standard fact [24J that the leading coefficient of the normalized polynomial 
Hk is equal to and we agree := 0. ^ = ^(C;A;) solves the following RH problem, 

which is a slightly modified version as the RH problem for Y satisfying the conditions ()3.1|) -( f3^2l) . 
but now corresponding to the external field x 2 instead of nV{x). 



RH problem for ^ 

(a) * : C \ R C 2x2 is analytic 

(b) For x G R, 



1/ ' 



(4.6) 



(c) ^ behaves as follows as £ — ► oo, 

* (C)= ( I+ ^(-j|g 7")+^(|) 2 )^^. asC-oo. (4.7) 

Note that the matrix on the right hand side of (|4.1|) looks similar to for the value of k = 0. 

4.2.2 Construction of the parametrix 

We will define the parametrix P of the following form, 

P(z) = E n , t (z)^(V^fnA^k)e- n ^ z ^(z - x^ t )~^\ (4.8) 

Here the parameter k in the model RH problem is, as before, equal to the non-negative integer 
that lies closest to v. The analytic functions E n ^ and f Ht t are still to be determined at this 
point. 
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Jump condition for P 

Since we assume E n ^ analytic in U x *, it does not effect the jumps of P. If f n t is a real conformal 
mapping in U x * with f n ,t(Xn t) = and f n t (x* n t ) > 0, it follows from the jump relation (|4.6p for 
that P satisfies the jump conditions 

P+(x) = P-(x) ( " 1 J , for x* n>t <x< < it + 5, 

/ p -2-Kiv |„ _ _* |2^„2n0 nit (z)\ 



This shows already that our parametrix P satisfies the jump condition we required in Section 
13.4,31 for any choice of E n ^ and f n j. The freedom we retain to define E n ^ and f n j, will be 
necessary to create a suitable matching of P with p(°°) at cJf/j;* . 

Matching condition for P 

In the double scaling limit where we let n — ► oo and t — > 1 in such a way that |t — 1| < M^p 1 , 
we would like the following matching condition to hold, 

P{z)P {po \z)~ 1 -> /, for z 6 dP x *. (4.9) 

The value of our second scaling parameter s = v will determine how good the matching between 
P and p(°°) can be. If v is an integer, we will have a good matching up to order 0(n~ 1 / 2 ). When 
v moves further away from an integer, the matching will become worse, but still reasonably good. 
Only when v is close to a half integer, the matching is not good anymore, and then we will only 
have matching up to order 0(1). 

Let us first define f n t similarly as in the case where t < 1 by 

f n<t (z) 2 = 2cf> n>t (x* n>t ) - 2<i> n!t (z). (4.10) 

Again we know by the definition of x* n t that Q n>t has a double zero at x* nt , or equivalently, that 
h n t has a simple zero at x* t . This implies using the definition (|3. 13|) of <p n t that the right hand 
side in (|4.10p has a double zero at x* n t . Consequently this again defines f n> t in a conformal way 
near x* , with 

fn,t(%n,t) = 0; fn,t( x n,t) = T^jyJ Q'n,t{ x *n,t) l ^ A ■ (4-H) 

As n — > oo and t — ► 1, it is clear from the discussion in Section [2] that 
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fn,t«t) = ^Qv^r + 0(t - 1). (4.12) 
Note also that the definition of f n< t ensures the existence of a constant C > such that 

Sn\f n ,t(z)\ > CV^, for z £ dU x *, (4.13) 



under the condition that we have chosen U x * sufficiently small but fixed. This means that, as 
n — > oo, we can use the asymptotic condition (14. 7ft when evaluating \f at \^nf n ,t(z). By (|4.8p 
we have that 

P(z) = P n X^)(/ + 0(n- 1 / 2 ))(v^/nX^)) fcCT3 e^ n ' t( < t)CT3 (^-< t )~ I ' CT3 J for z E dU x *. 

(4.14) 
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This behavior suggests how we should choose the analytic pre- factor E n f. If we take E n ^ as 
follows, 

E n , t (z) = P^\z)(z - x^ t r 3 e n ^<^(^f n , t (z))- k ^, (4.15) 

one checks directly using the definition (|3,29p of p(°°) that E n j is analytic in U x * \ {x^ t } with 
a removable singularity at x* t . 

Inserting this definition of E n t into (|4. gives us the following behavior of P for z G dU x * 
in the double scaling limit, 



P(z) = P^\z)(z - x*)Y^e n ^<>H^f n , t {z)y 



ka 3 



ikl 



x (^WJ^e^tO^-xy-"". (4.16) 
We can express (|4.16p in the following more convenient form, 

P(z)P°°(z)- 1 = (7 + ©(n- 1 / 2 )) ^^(z)- 1 (4.17) 



If E n t would be bounded uniformly in n and t for z £ U x * , this would provide a good matching. 
However, for general values of v this is not the case. We can write 

E n , t (z) = E n .t{z)e n ^<>^- ka \ (4.18) 

where E U; t is bounded uniformly in n and t, since it depends on n and t only through the points 
o«n,ti b Ut t, and x* t , which vary smoothly with n and t. The quality of the matching now depends 

on the asymptotic behavior of e a ^ n ' t ' yXn ^ a ' i 'yfn ka ' A . We deal with this behavior in the following 
proposition. 

Proposition 4.1 In the double scaling limit where n — ► oo and t — » oo in swc/i a way t/iat 
It - II < with 

n f 1 * 1 

i/ = 2(t-l)- / = <foeR+, (4.19) 

we have that 

Proof. We show that 

n<t>n,t( x *n,t) = logn + 0(1), 

from which the proposition follows directly. Checking the variational equality for the measure 
tpi n> t, starting from (jl.6p . learns us that t// n ,t is the equilibrium measure in external field V, 
with mass t(l — m n> t). Now it follows from a formula by Buyarov and Rakhmanov [3] that 

/t(l— m„ t t) 
u n ,sds, (4.20) 

where uj n j denotes the equilibrium measure minimizing the unweighted logarithmic energy 
= JJ log | x \ — ] dw{x)dw{y) 
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among all probability measures supported on the support [a' n t ,b' nt ] of the weighted equilibrium 
measure [i n ^t- In other words the equilibrium measure in an external field can be realized as an 
integral of unweighted equilibrium measures. Further it is known [24] that 

dcj n ,t(x) = 1 dx. 

TTyJ{x - a' n t )(b' n t - x) 

On the level of densities, (|4.20p means that 

th nit (z)yj (z - a' ntt )(z - b' n t ) - hn tl (z)y/(z-a)(z-b) 



/t(l-m„,t) i 
, =ds- (4.21) 

yj ( Z ~ a 'n,s)( z - K,s) 



Integrating this identity from b' nt to x* nt gives us, using the definition (|3.13p of 4>n,t, that 



b't 



t<f> n ,t( x n t t) - 0n,l(a?n,t) + / h nA (z)y/(z - a)(z - b)dz 



dsdz. (4.22) 



Using the fact that (j> n ,i(x*) = 4>' n ±{x*) = and estimating the third term on the left using the 
behavior of the special points 

x* ntt = x* + 0(t-l), b' njt = b + 0(t-l), 

we obtain that 

t<t>nMn,t)= / , ^=dsdx + 0(t-lf/ 2 , 

Kt Jl ^{x-a' n>s ){x-V n>s ) 
Now using also the facts that ra n ^ = ^ = C(n _1 ) and that a' nt = a + 0(t — 1), we find that 

n , t « it ) = (t-l) f 1 dx + 0{n- 1 ), 

Jb \/{x — a){x — b) 

and consequently by (j4.19H we find that 

n<t>n,t( x n,t) = o logn + 0(l), 
which proves the proposition. □ 
It follows from the proposition that, in view of (|4.18p . 

E n<t {z) = C(y/n lAl ), with A = v - k, 
and this gives us by (|4.17p a matching as follows, 

P(z)P 00 (z)~ 1 = 1 + 0{n~ l / 2+ \ A \). (4.23) 

This ends the construction of the local parametrix. 

Indeed we see that the quality of the matching is good when v is close to an integer (or 
when A is small) , but is getting worse when v approaches a half integer (or when A approaches 
±1/2). However, even near the half integers, the parametrix fits sufficiently well to obtain the 
asymptotics (|3.40p for R, which is necessary to prove Theorem 11.61 
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5 Universality of the eigenvalue correlation kernel 



In the previous sections we found asymptotic for R in a suitable double scaling limit where 
n — ► oo and t — ► 1. Reversing the transformation S i— ► P and using the explicit formula (|3.2ip 
for the kernel in terms of 5, will enable us to find asymptotics for the eigenvalue correlation 
kernel and to prove Theorem 11.21 

For x, y £ M D U x *, we have by (|3.36[) that R(z) = S{z)P(z)~ l , with P the local parametrix 
constructed in Section [H Inserting this into (|3.2ip gives us the following identity, 

(x - x* f ) u ,(y - xl f )le n ^+^e n ^+^ , . (\\ 

KnMv) = - 2m(x - y) ( ° ^ P + WW (qJ " C 5 " 1 ) 

We will now use the asymptotics (|3.40p for R and the explicit formulas for P in order to find 
asymptotics for K n< t. 

5 . 1 Proof of Theorem Q] 

We start with the proof of Theorem 11.21 in the one-interval case where t < 1 and i/ = 0. We 
prove that the re-scaled eigenvalue correlation kernel is trivial in this case. 

Proof of Theorem 11.21 if t < 1. We first recall, from formula (|4.ip . that the structure of the 
local parametrix is as follows, 



P(z) = P(°°\z) (I *V for zeU x *, 



where * denotes a matrix-entry for which the precise value is unimportant. Now plugging this 
parametrix into (|5.1|) leads to the following equation, 

= 2vn(x - y) ( ° l) P^HyT'mR-H^P^W (q) ■ (5-2) 

Because of the uniform asymptotics (|3.40p for R and the analyticity of R in the disk U x * , we 
have that 

R~ 1 (y)R(x) = I + O C^jf) as x,y ^ x* and rw oo. (5.3) 

Since p(°°) is analytic in U x * (if v = 0) and uniformly bounded in n and t, we obtain for n — > oo 
and x,y ^ x* , 

e n<t> n ,t,+ (x) e n<t> n ,t,+ (y) / /i\ f T — v\\ 

- ( (o ^""w'^K^ ^)) 

P "</'n,t,+ (a;) (= n.0n.t,+ (j/) 

Now recall from (|3.7p (where in addition we should note that, since //t = pt as t < 1, inequality 
holds near x* as well) and (|3.15p that <f) n t( x ) < for x > b' n t as t < 1. This is already sufficient 
to conclude that K n t {x,y) is bounded, which clearly implies, for 

X = X* + -T7TT, X = X* + 



(cn)V2' (cn) 1 ^' 
that 

—^—K n (x, y) = 0(n- 1 ' 2 ) as rw oo. (5.5) 
This proves the theorem in the case t < 1. □ 
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Remark 5.1 It is worth noting that we did not use the fact that we re-scaled with a factor 
(era)" 1 / 2 . If we would re-scale by putting x,y = x* + 0(n~ 7 ), for any 7 > we obtain directly 
from d53J that ±K n>t (x,y) = 0(n~ 7 ). The only reason why we did choose this particular 
scaling, is that it is the 'right' scaling for t > 1. 

For t > 1, we have that v > 0, which leads to a local parametrix which is somewhat more 
complicated. This makes the proof of the theorem slightly more involved. 

Proof of Theorem 11.21 if t > 1. We recall from Section [4] and (|4.8h in particular that the 
local parametrix P has the form 

P(z) = J E n , t (z)^(^/ n , 4 (z);A ; )e-^ t ( z ) CT3 (z-< it )-^. (5.6) 

Inserting this formula into (15. ip leads us to the following identity, 

K n ,t(x, y ) = 27Ti( l_ y) (0 1) ^{V^uM^Kl^Riy) 



x Br\x)E njt {x)*+{yfcU, t {x)\k) r\ . (5.7) 



As in the case where t < 1, we derive from the uniform asymptotics (|3.4U|) for R and the 
analyticity of R that the following holds in the double scaling limit, 

R-\y)R(x) =1 + ( J /2 ~ f A| ) , far a, y-z*. 

The structure (|4.18p of E then implies together with Proposition 14.11 that 

E nyt (z) = 0(V^ Aaa )E n , t (z), 
and this yields 

E-\{y)R{y)R- 1 (x)E nyt (x) = I + O ( (x - y)n) A \) , iovx,y^x*. (5.8) 



Let us now re-scale the variables x and y by putting 

X = x * + T^m> y = x * + T^> withc = f[(x*) 2 = ^=q'v(x*) 1/2 , 
(en) 1 ' 2 (cn) 1 ^ V2 

where u,v £ [-M, M] for some arbitrary large constant M > 0. Since we know by (|4.12p and 
from Section [2] that 

f' nt (x* nt ) 2 = c + 0{t — l), x* n t = x* + 0(t — 1), as n — > 00 and t — > 1, 

we have that 

Vnf n ,t( x ) = u + 0(Vn(t - 1)), VnfnAv) = v + 0(Vn(t - 1)). 
Using the above estimates, (I5.7P reduces to 
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If A remains a fixed distance away from ±1/2, it is now a straightforward calculation using the 
definition of ^ to check that 

lim - 1 1/0 K„ t ( x* + - M 1 /0 , x* + 



{en) 1 / 2 rhl V (cn) 1 /2' ( m )l/2 

fc e - " 2 +" 2 H k (u)H k __i(v) - H k (v)H k _i(u) 
2 u — v ' 

which completes the proof of Theorem 11.21 □ 

5.2 Proof of Theorem 11.61 

In order to prove Theorem 11.61 we need to find more accurate asymptotics for the eigenvalue 
correlation kernel. We will do this by exploring in more detail formula (|5.7p . For the proof of 
Theorem II .21 it was sufficient to approximate E~\(y)R(y)RT l (x)E n . t {x) by the identity matrix. 
To arrive at Theorem 11.61 we need to be a little bit more careful, and we need some better 
estimate, compared to (|5.8p . From the asymptotic formula for R we obtained in Section [3.51 we 
recall that 

R{z) = / + 0(n- 1 / 2+ l A l), (5.11) 

uniformly for z € C \ E_r, in the double scaling limit we considered. 

Using (|4.18p and Proposition 14.11 we can refine (|5.8p in the following way as x, y — > x* in 
the double scaling limit, 

E^ t (y)R(y)R- 1 (x)E njt (x) 

(x - y)n A + 0(x - y), as A > 0, 



= < 











/ + 

< 






:) 






fo 




I — 






:) 



(5.12) 



(x - y)n~ A + Oix - y), as A < 0, 



where c nt are some sequences that are bounded in the double scaling limit. 

Proof of Theorem 11.61 We start by picking up the exact formula (|5.7p from the proof of 
Theorem 11.21 We now follow the same calculations as in the proof of Theorem 11.21 but with 
(|5.12p inserted into (|5.7p instead of the less accurate approximation we used before. For A > 
and with 

it is straightforward to check that (|5.9|) should now be replaced by 

K " AX - V] = 2xi(t-y) ( ( ° (J) 

+ 2* <t nHr.-y){0 °) * + («; k) Q + O ( ( " -J»;° g " )) . 



We have already computed the first term on the right hand side in the proof of Theorem 11.21 
leading to the Hermite kernel K GUE (ii, v; k). Using (|4.5p . we can easily compute the second term 
as well, 

(o i) *:?(«;*)(; o) (J) = ^e-^H k (u)H k (v), 
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which gives us the formula 

The Christoffel-Darboux formula for orthogonal polynomials allows us to rewrite the right hand 
side of this equation as 

e-^E WW + <,^-^H*(.>**«ji^ + O (^) , (5.13) 



j=0 

which is equal to 



(1 - A+ t )K^ u >, v; k) + \t t K^(u, v; k + 1) + O ( ) , (5.14) 



log n 



with 



+ _ + y^fc! 1 

A n,t ~ C n,t 2 k n l/2-A' 

For A < 0, we obtain in exactly the same way that 
1 



.^^K n , t (x,y)=K^(u,v;k) 

Again using the Christoffel-Darboux formula, this can be rewritten as 

—L^tfo y) = (1 - A" ,)IfC GUE ( U , «; fe) + A" ,K GUE (u, «; k - 1) + , (5.17) 

with 

2 fc+1 1 

K,t = c n^ (A; - l)! n l/2+A- ( 5 ' 18 ) 

This proves (11.14h . and the estimates (j!.15[) and (|1.16j) for follow from f)5. 15|) and (|5.18j) . 
This completes the proof of Theorem 11.61 □ 

Remark 5.2 One can obtain explicit formulas for the constants c nt by computing the sub- 
leading term in the asymptotic expansion for R. This would provide explicit formulas for X^ t , 
which are however rather complicated. Since we have not been able to simplify those expressions 
considerably, we feel it is not very useful to give the involved computations leading to those 
formulas here. 
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